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Abstract 

The result of performing integrations over connection type variables in the 
path integral for the discrete field theory may be poorly defined in the case of 
non-compact gauge group with the Haar measure exponentially growing in some 
directions. This point is studied in the case of the discrete form of the first order 
formulation of the Einstein gravity theory. Here the result of interest can be 
defined as generalized function (of the rest of variables of the type of tetrad or 
elementary areas) i. e. a functional on a set of probe functions. To define this 
functional, we calculate its values on the products of components of the area 
tensors, the so-called moments. The resulting distribution (in fact, probability 
distribution) has singular ((^-function-like) part with support in the nonphysical 
region of the complex plane of area tensors and regular part (usual function) 
which decays exponentially at large areas. As we discuss, this also provides 
suppression of large edge lengths which is important for internal consistency, 
if one asks whether gravity on short distances can be discrete. Some another 
features of the obtained probability distribution including occurrence of the local 
maxima at a number of the approximately equidistant values of area are also 
considered. 
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1 Introduction 

Suppose we are studying discrete version of a theory with (one of) the field variable 
being of the connection type. Example is Wilson formulation of the QCD[1]. The 
discrete version of the connection is not an element of the Lcc algebra but an element 
of the gauge group itself. In the path integral when integrating over the discrete 
connection type variable one is faced with integration over Haar measure on this group. 
In the case of (the first order formulation of) the Einstein gravity the group is that of 
the local rotational symmetry of the spacetime itself, i. e. S0(3,l). This is noncompact 
group while the Haar measure exponentially grows with Lorentz boost angles which 
together with Euclidean angles parameterize the discrete S0(3,l) connection. The 
corresponding integrations in the path integral should be made with care. The result 
of performing these generally does not exist as usual function (of the remaining variables 
of the tetrad type) but it can be given sense of the generalized function, or distribution 
(in fact, probability distribution for the remaining tetrad type variables). In fact, 
the latter turns out to include (5-f unctions and their derivatives, but it also have a 
regular part being usual function. Such structure of this distribution does not allow to 
continue it to the Euclidean region (5-functions do not exist for complex arguments). 
This correlates with the fact that Euclidean path integral for pure Einstein gravity 
requires careful definition because of the unboundedness of the action from below. 

Both discrete and continuum versions of the first order Einstein gravity are natu- 
rally formulated using decomposition of S0(3,l) element into two mutually complex 
conjugated elements of SU(2,C) or S0(3,C) (complex- valued) . For definiteness, one 
might imply the following notations and sign conventions concerning splitting the ten- 
sors into (anti-) self dual parts in the Minkowsky spacetime. Suppose there is S0(3,l) 
matrix (being connection type discrete variable), 

n = e^v{^''El, + ,p'Ll,). (1) 

Its generator is expanded over the set of independent generators, 

Ekah — —^kab, L^ab = 9ka90b " 90a9kb {9ab = diag(— 1, 1, 1, 1), £123 = +!)• (2) 



We denote 



^fcab — —^kab ± i{9ak90b — 9a09kb) 



(3) 



so that 
then 

^ = exp(^^^^^^,). (5) 

In the continuum hmit the generator of Vl corresponds to ujf'dx'^ where uj'^ = — cu^" 
is the infinitesimal connection which is thus expanded into self- and antiselfdual parts 
additively, 

For a triangle spanned by the two 4- vectors II, we can define bivectorv"'^ = |e"''cd^i^2- 
This variable (or antisymmetrized tetrad bilinears in the continuum theory) or used 
somewhere in the present paper its analytical continuation to arbitrary antisymmetric 
area tensor split in the same way, 

^ab ^ +^ab ^ -^ab^ ^ab ^ l^ab ^ (7) 

In particular, 

2^ o ^ — V o V ± iv * V. (8) 

Here AoB = \AabB''\ A*B = leabcdA"^ B"'^ for the two matrices A, B. The ^-parts 
map into three-dimensional vectors 

^Vab = ^^'' ^kab, 2 = -eklmV^"" ± i{VkO - Vok)- (9) 

For a bivector 2 % = ±ili x I2 — Z1/2 + hli- Additional overall i is here due to the fact 
that v"''' is dual area tensor. Besides that, 

V = 2^o^. (10) 

The is (—1) times the square of the (real for the spacelike triangle) area. 

In the case of the continuum Einstein gravity the most general first order form of 
the action is that by Hoist [2, 3], 
1 r 2 

-^Hoist = g y (ea6cde^ej;+ -exce^,d)^^^''''{d^^p - d^uj^+ u^Up- ujpUJ^y^d'^x, (11) 

which generahzes the Cartan-Weyl form of the Einstein action and reduces to | / 
y/—gd'^x in terms of metric gxf^ — e\ea^ if the infinitesimal connection ujf' is expressed 



in terms of the tetrad e" via equations of motion for a;"^. Upon splitting cu into (anti- 
)selfdual parts the S'noist sphts too and is easily seen to have the form (1 + 2/7) '^Scw + 
(1 — i/'y)~Sc-w where ^^cw arc ( ant i-) self dual parts of the Cartan-Weyl continuum 
action (that is, of (ShoU at 7"^ = 0), 7 is known as Barbero-Immirzi parameter [4, 5]. 

The discrete minisuperspace formulation of general relativity on the piecewise fiat 
manifolds or simplicial complexes is known as Regge calculus [6, 7]. Invoking the notion 
of discrete tetrad and connection first considered in Ref. [8] we have suggested in Ref. 
[9] representation of the Einstein action | / R^—gd'^x on this manifold in terms of area 
tensors and finite rotation S0(4) (S0(3,l) in the Minkowsky case) matrices, and also 
in terms of (anti-)selfdual parts of finite rotation matrices. For the latter we write 

= !L V arcsm == . (12) 

Here are area vectors (9) of the triangle cr^, in the Minkowsky case Q^-s is con- 
nection SO (3,1) matrix on the tetrahedron which we call simply connection, R„2 
is curvature matrix on the triangle a^, holonomy of fi's (some product of f2's). For a 
3- vector v and a 3 x 3 matrix R we have denoted v*R= ^v°'R'^'^€abc, and for "^R^^ , the 
(anti-) self dual part of R^2, we have used adjoint, S0(3) representation (to be precise, 
S0(3,C) matrix). 

These representations result in the same | / Ry^^d'^x on the piecewise fiat mani- 
fold in terms of the purely edge lengths upon excluding rotation matrices by classical 
equations of motion (that is, on-shell). Taking into account that in the Minkowsky 
case +5" = ( ~S)* we can write out the most general combination of "^S", ~S which i) is 
real and ii) reduces to | / Ry/^d'^x in terms of the purely edge lengths on-shell, as S 
^ C+S + C*-S where C + C* = 1, that is C = [1 + i ■ (real parameter)]/2. We see 
that there is direct analogy with the discrete case if we write C = {1 + i/'~f)/2 where 
the discrete analog of 7 is denoted by the same letter. We assume < 7 < 00. Thus 

Let us write out a discretized functional integral with this action, / exp{iS)Dq, 
q are field variables (some factors of the type of Jacobians could also be present). 
Functional integral approach in Regge calculus was earher developed, see, e. g., Refs. 
[8, 10, 11]. Suppose we have performed integration over rotation matrices and are 
interested in the dependence of the intermediate result on area tensors. Of course. 
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different (components of) area tensors are not independent, but nothing prevent us 
from studying analytical properties in the extended region of varying these area tensors 
as if these were independent variables. Namely, consider integral 



Matrices can be parameterized by complex vector angles ^(f) — (p^itjj (rotation 



A special feature of the integrals like (14) is that these arc generally not usual 
functions oi ^v, ~v — {~^v)* but rather generalized functions or distributions. This 
is because the Haar measure X>Q grows exponentially with imaginary angles (i. e. on 
Lorentz boosts), and integral over it may diverge as usual function. This is consistent 
with the sense of the result of intermediate integrations in the path integrals over 
connections as some objects to be further integrated over area tensors. Therefore in 
what follows we just treat these integrals as distributions, i. e. study integrals of these 
expressions with some probe functions. For the probe functions taken as products of 
components of % we can give sense to these integrals (called moments in this case). 
This allows to define distributions of interest. Remarkable is that the singular (5- 
function-like) part of the latter turns out to have support outside the physical region 
Im v"^ — 0. So this distribution in the physical region is usual function. 

Before proceeding with analysis of such integrals it is useful to maximally separate 
out the gauge degrees of freedom contained in the connection variables ilg-s. The 
matrices R^2 possess more physical sense. Only part of these are independent; those 
on certain set of the triangles are functions of other i?o-2S, a"^ ^ T due to the Bianchi 
identities[6]. A natural step is to reduce integration element Ho-s ^^^^o-s to IIo-^^jf^-Ro-^ 
times the product of PQo-s's on a certain set Q of cr^s which absorb the rest of the 
gauge degrees of freedom, W„z^qV^„z. 

Let us ilhistratc this by explicitly expressing O^-^^s in terms of {Ra^ '■ o"^ ^ J^} and 
{f2(^3 : CT'^ G ^} for the simplicial structure regular in certain direction. This structure 
specifies the idea of constructing the piecewise flat 4-geometry of separate piecewise 
flat 3-geometries[12]. The simplest such structure is widely used periodic simphcial 
decomposition when the spacetime is divided into 4-cubes and each 4-cube is divided 




(14) 



by the angle J around the unit vector ^4>l \l =^^). 
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into 24 4-simplices[13] by diagonals emitted from certain vertex of the cube. This 
structure is uniform w. r. t. the choice of vertex and similar along any of the four 
directions along the 4-cube axes. (If it is important to have finite number of vertices in 
lattice analysis, one might easily pass to the torus topology by imposing cyclic boundary 
conditions along corresponding directions.) We use our general notations[14] suitable 
for this structure as a particular case, see fig.l. 



Namely, we consider the sequence oi t — const 3D simplicial 3-geometries, call 
these leaves, separated by, say, the change At of the coordinate t such that for each 
vertex i in any of the leaves there arc both its image and pre-image in the two 
neighboring leaves a,t t — At and a,t t + At. The coordinate t is denoted as the (world) 
time although it well can be viewed as space coordinate. The points of the leaf at t are 
denoted by the letters i, k, I, . . .. Let us denote a simplex by enumerating its vertices 
in round brackets. Besides that, we assume the same scheme of connection of different 
vertices by links in the leaves. Let us call these links and other simplices like, e. g., 
tetrahedron [iklm] completely contained in the leaf, the leaf simplices. The full 4D 
complex is supposed to result from these 3D leaves by triangulating the space between 
each two neighboring ones by new, diagonal and t-like links. The link (ii"*") or 
which connects a vertex with its image in the neighboring leaf, will be called t-like link. 
Any other simplex like, e. g., the tetrahedron (ikk^l) containing such link, will be 
called t-like simplex. The diagonal simplices are the remaining ones, e. g. tetrahedron 
(ik'^lm). 

The spacetime is thus divided into the "i-hke" four-dimensional prisms. A face of 
the latter, the three-dimensional prism is just shown in fig.l. 

Then R„2 on the leaf and diagonal triangles could be taken as independent ones, 
and i?o.2 on the t-like triangles are functions of them via Bianchi identities. The Qg-s 
on the leaf and diagonal tetrahedrons could be taken as purely gauge ones. (All 
this in particular implies that the number of the leaf and diagonal triangles matches 




Figure 1: Fragment of the 3-prism. 
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the number of the t-\ike tetrahedrons. This is indeed so.) Indeed, consider a chain 
of expressions for the R„2 on the leaf and diagonal triangles sequently filling in a 
3- prism, see fig.l. 



R{ikl) — ■ 


■ ■ ^^{i-ikl) ■ ■ 


■ ^{ik+kl) ■ ■ ■ 




R{ik+i) = ■ 


■ ■^''{ik+kl) ■ ■ 


■ ^(ik+l+l) ■ ■ ■ 


(15) 


R{ik+i+) = ■ 




■ ■ ^{i+ik+l+) ■ ■ ■ 





The dots in the expressions for R mean matrices on the leaf and diagonal tetra- 
hedrons. From Eqs. (15) using some initial values for fls on i-like tetrahedrons at 
t = -oo we can step-by-step express . . . ^^^+^^+1+) ^{ik+i+i) ^{ik+ki) ^{i-iki) 
. . . (the arrow means " in terms of ) as functions of Rs and Qs on the leaf and diagonal 
triangles. In particular, knowing these fls we can find the curvatures on t-like triangles, 
the products of these fls, e. g. 

^(i+ik) - ^I'ii+ikU) • • • ^ ^{i+ikh) ■ J 

(This is just what shows up inside the 3D section as a curvature on the link {ik).) Here 
€-{iki)m = il is some sign function. Besides that, we can step-by-step rewrite integration 
element in another variables passing from the Vts on t-like tetrahedrons to Rs on the 
leaf and diagonal triangles as ■ ■ .^{ik+ki) R{iki), ^{ik+i+i) R{ik+i), ^{i+ik+i+) 

R(ik+i+), On each step we have typical relation between certain matrices R and 

Q as = rii?r2 where SO (3,1) matrices Fi and Fi are products of another matrices 
and R temporarily treated as constants. Then VQ — VR due to the left-right 
invariance of the Haar measure. We thus arrive at 

\[vn,s = n VR,2 n vn,^ (17) 

where T is the set of i-hke triangles and Q is the set of the leaf and diagonal tetrahe- 
drons. Thus integration in (14) over connections could be replaced modulo integration 
over some connections taken as "gauge" ones by that one over independent curvatures. 

The present paper simply calculates the simplest integral of the type of Eq. (14). 
Outside of nonphysical singularities this function exponentially decays at large areas. 
Some features of this function (Euclidean version, occurrence of a number of local 
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maxima), consequences for the general case, for the measure on the edge lengths, are 
discussed. 



2 The moments of probability distribution 

2.1 Defining the moments 

As mentioned above, we consider Eq. (14) as function of arbitrary ^v, "v — {'^v)*. 
Consider moments of A/" (integral with powers of %). As also mentioned, Af contains 5- 
functions which cannot be continued to/from the complex arguments. However, in the 
definition of moments the contours of integrations can be deformed to complex plane. 
Namely, integration contour over i/) considered as complex variable is deformed from 
real to imaginary values. The become independent real variables, and integration 
over Vfl splits, 

vn^v-^v-n, M^+N-N, p^ = ^^^^t^d3±0. (18) 

We define M in the region where I — (l±i/7)%/2 are real and then continue to % 
of interest. We are faced with integrals of the form 

jv^R...j e*'^("^)r...rd^Z (19) 

where I — VP, n — l/l, g{x) is odd analytical in the neighborhood of x = function, 
g{—x) — —g{x), such as principal value of arcsinx of interest or simply x. The vector 

r"' = e^i,c^R^'^/2 = {^(j)°'sm =^0)/=^(/). Existence of integrals (19) easily follows at g{x) = 
x: we simply get derivatives of 5-functions 5^{r) which are then integrated over 

Finiteness is provided by analyticity of this measure at r ^ 0, P^i? = (cq + cir^ + 
C2{r'^)'^ + ...)d^r". Dots in JVR... in (19) mean possible dependence on R of factors 
provided by Rc^2 on other triangles due to the Bianchi identities. 

For the more general case g{x) ^ x (19) is also finite. Again, consideration may go 
through appearance of 5-functions at an intermediate stage. Namely, special structure 
of exponential (19) allows to extend integration over I to the whole range (— cxo, +oo). 
This is only possible because formal substituting / — > — / is equivalent to n — > — n due 



to the oddness of g{x). This results in 5-functions of g, 

+00 

—00 

= ^{2n){-iy+^ J S^^+^\g{nr))n''K..n''^d^n. (21) 
Apply S'^'^'^\g{x)) to probe functions, 

da; 



/ 



5^^-''\g{x))f{x)dx 



d \^'+^ 



(22) 

5=0 



we find it be expanded over 5^^+^\x), S^^\x), ... , S^^^'^°'^'^^+^\x), 5^^^'"°'^^^\x). For 
each term in this expansion in the RHS of the formula (21) we can use analogous to 
Eq. (21) relation (at g{x) — x) but read from right to left for the backward conver- 
sion to the expression having form of the Fourier transform of power function of I, 
I. The leading term 5(-'+^^(nr) in this expansion is converted to / exp(iZr)Z"^..Z"^d^Z 
= {27cf{-iyda,...da^6^{r) (the LHS of the Eq. (21) at g{x) = x). This being inte- 
grated over V'^R is finite. Subsequent terms 5^^~'^''\nr), < k < [j/2] ([j/2] is integer 
part) are converted to (— i)-'"'"^(?oi...c?a^27r^(r^-|-£^)*^~^/^[(2A;)!]~^. Here intermediate reg- 
ularization £ ^ of the singularity at r = is convenient to detect possible formation 
of the contact terms via 9^(l/r) = — 47r5^(r). Finiteness is provided by analyticity of 
the measure at r = 0, 

m=0 

9{y) is Heaviside step function. Convergence of integral with /„(r)^(l— r), n — [j/2]—k, 
is seen immediately, and that of integral with P„(r)^(l — r) follows upon multiple 
applying integration by parts. 

2.2 The moments of the simplest integral over connections 

Above we have shown possibility to define the (nonabsolutely convergent) expressions 
for the moments in a finite way. If simplified integral is considered. 



jVo = y exp ^ 



i\ /—t; . ^v*~^R / i\ / — - . V* R 
H-— V "t;^ arcsm — . -|- 1 v "i;^ arcsm — . — 



2 



VR 
(24) 



(a dependence on R of R^^ on other triangles due to the Bianchi identities is not 
taken into account), calculation of arbitrary moment of it can be performed in closed 
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form. It is sufficient to integrate in (19) and then over T> the scalars P'^. Substituting 

f{x) = /(2™)(0)a;2™/(2m)! and j = 2k into (22) gives {d/dgf''+^[x^"'+\g)]/{2m + 
l)!|g=o for the coefficients in the expansion of 6'''^''^'^\g{x)) over 6^'^"^\x)s. Using (21) 
read from right to left we reduce / 5^'^"^\nr)d'^n to / ex.p{ilr)P"^~'^d^l. Summation 
over m with the coefficients found allows to find / exp[ilg{nr)]l'^'^d^l so that 



'A' 

A9. 



2A;+3 



E(-ir 



.m=0 



(2m + 1!) 



5=0 



2fc+2 



with the "generating function" 



dx 



(25) 



J 9=0 



lix) 



d^i 



/ 



jir cos(a^0 

/2 



dH 



.r/(-^-l)^ = .ln(l + Vr^). 



(26) 



We have extended summation to infinite number of powers of x keeping in mind that 
upon applying (d/dg')^'^"^^(-)g=o only finite number of terms are active. 



3 Restoring probability distribution from the mo- 
ments 

In the simplest case x = g the moment (25) is 
At X = sin g the moment is 

M [cos5ln(l + cos^)]^^o. (28) 

It is not difficult to find out density of distribution giving these values on monomials 
(perform a kind of Mellin transform). Appropriate entries in the table of integrals are 

1 2 

g = - / chglKo{l)dl, (29) 
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Ko is modified Bessel function, and 

1 1 



|sinc/-- + -cos^ln[2(l + cosc/)] = y j^-^ 



Thus we find at g{x) = x 



I chgl 
+ 1 shnl 



dl. 



(30) 



dH 



(31) 



where 



Ki^H) = [ e-'-h,^ = / exp — dip 

J chrj J \ sin (/? y 



(32) 

\ Bin 0? / 

'0 \ r / 

is integral of Ko{l) — —[Kii{l)]' , Kii{oo) — 0. This had appeared in our work[15] in 
3D SO (3) Regge calculus. 

At g{x) = arcsinx of present interest 



+ (21n2-4)/(l)-4/'(l). (33) 



shTrZ P + 1 7r2 

The RHS in both cases is normalized so that appearing there integral be 1 at / = 
1. For the integral (24) we deform integration contour and pass from Z to = +v = 
2/(1 + i/7)-i and v*, 

I AAo(v, v*)f{v')h{vyd'vd'v* = i^{f)m\ (34) 

where v = \/v^ (v^ + v*^ = — Tr(t'^)) and 

z f (l/7-^)v/2 /(v^)dS 
^^■'^ 2 J (l/7-?)2vV4 + lsh[7r(l/7-i)v/2] 
+47r(l + i/7)-'[(ln 2 - 2)/(4(l + i/7)-^) - 8(1 + i/7)- V'(4(l + ih)-% (35) 

Here additional to integral terms have support at nonphysical points = 4(1 +i/7)~^. 
Thus in any region excluding these points (e. g., in physical region Imv^ = 0) we have 

2 



A^(v,v*) = 



1 



Hi-)- 



(36) 



Looking at the logarithmic expression for generating function (26), it is interesting 
to ask to what extent our definition based on continuation from Euclidean-like region 
relies on the compactness of SO (4). Alternative way to define the moments of (24) 
might proceed via deforming integration contours so that (p be imaginary. The ^0 
become independent imaginary variables (pure Lorentz boost angles) , and integration 
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over VR splits. The r is imaginary, we define moments on imaginary I and find I{x) at 
imaginary x. Logarithmic divergence at r — oo appears. Let the cut off for r in (26) 
after replacing r ^ ir be vq. Then continue I{x) to real x. The I{x) (26) is modified 
by adding constant to logarithm, 



7(x)^7rln|ii^^^2|. (37) 



1 + Vl + ^o, 

However, it is not difficult to see that this constant only modifies the coefficient at 
/(4(1 + i/7)~^) in (35) and thus does not affect the result (36). 

In physical spacelike region v^ = v*^ = — |vp, |v| is module of the triangle area, the 
Eq. (36) behaves as cxp(— 7r|v|) at large |v|. In physical timelike region v^ = v*^ = |vp, 
this behaves as exp(— 7r|v|/7) at large |v|. For comparison, replacing arcsinx by x as 
in (31) gives exp(— |v|) and exp(— |v|/7), respectively. 

An alternative way to define moments does not include transformation to the Eu- 
clidean signature [16]. This approach encounters more singularities and requires care- 
ful mathematical analysis. The resulting probability distribution coincides with the 
present one. 



4 Discussion 



4.1 Qualitative explanation of exponential suppression 

The above exponential suppression over areas arises already in the model integral 

J giv/^shv-^^ _ 2Ko{V^). (38) 

— oo 

Here \/— v^ is modeling module of the spacelike area, ip is modeling Lorentz boost 
angle. This behaves as exp(— V— v^) at large v^. Nonzero 7"-'^ mixes spacelike and 
timelike area components and leads to exponential suppression also in the timelike 
region. Taking integrals over connections (14) also reminds calculating Fourier trans- 
form of 1, although on curvy (group) manifold. While usual Fourier transform gives 
5-function, limiting case of exremcly rapidly decreasing exponent, calculations on curvy 
manifold lead to broadening this (5-function. This is quahtative "physical" explanation 
of the exponential suppression for areas. In particular, when g{x) — arcsinx, rotation 
angles enter exponential in the form close to the angles themselves. Then connection 
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integral is more close to the usual Fourier transform and provides larger exponential 

suppression for areas than the integral at g{x) = x (when exponential contains angles 
in the more nonlinear manner, as hyperbolic or trigonometric functions). 



4.2 Euclidean version 

Our definition of the integration over connection suggests that the Feynman path in- 
tegral itself is defined (on area tensor monomials, in a theory where these are treated 
as independent ones). However, direct evaluation of the Euclidean path integral fails 
because of the unboundedness of the action, although there are indications that the 
problem can be handled in the lattice theory, via appearance of a non-trivial ultra- 
violet fixed point (see, e. g., review [17]). We could modify definition. Namely, the 
Euclidean version of the path integral with exp(— jS'e), with the Euclidean signature of 
space-time and SO (4) connection could be defined if integration contours are deformed 
via formal change of the integration variables v —iv[18]. Namely, supplying (real) 
Euclidean variables with subscript E, follows formally from S (13) appearing in (14) 
by substitution ^0-2 ^o-^e, ^a^E, 7 ^7e (the latter because the coeffi- 

cients (1 ±i/7) should become real). Now 3- vectors '^v^^e and "v^^e are independent 
variables, as well as the elements of real SO (3) "^o-sg and "Hq-se are. The appropriate 
path integral can be defined on the area tensor monomials by deforming integration 
contours to purely imaginary tensors, ^o-ze ^a'^E so that the monotonic expo- 

nent exp(— S'e) becomes oscillating one. After finding N or, more exactly, Ae in this 
region, we should return to the original region of tensor values, ^o-^e ~^ i ^o-^e- Along 
the fines of sections 2, 3 thus defined Eucfidean version of Eq. (36) reads 

Moe{ +ve, "ve) = (39) 



i (i+l)^E + l Sh [I + i Sh [f (4-l)^E 



Here = y ^1- physical region {v^ is hivector) "K^e = ~ve = ve- Using here 
Minkowsky 7, 7e — > — ^7, we reproduce suppression exp(— ttve) of the triangle with 
area ve in the spacelike region, as it should be for the Euclidean version. To continue 
Euclidean result to the timelike region, one should know it analytically beyond the 
physical region +ve = ~ve. 
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4.3 General case 

In general, we have suppression of the triangles of large areas at least for part of the 
triangles, such that the resting part might be taken as those passing through the links 
being discrete analogs of lapse-shift. 

In more detail, when returning from the simplest integral (24) to the general one 
(14) the moments (25) get naturally generalized to the integrals of jV with monomials of 
the components for the whole set of the 4-simplices cr^. Important new point is that 
this integration should not be over alld^v^'^s. This is especially clearly seen in the case 
of the simplified action with the function g{x) = arcsinx replaced by g{x) = x. Then 
integrating exponent in (14) over (i^Vfj2 with the components v'^ for a given o"^ leads to 
the derivatives of the (5-function 5^(i?o.2 — 1(^2)- The set of holonomies {i?o-2 cr^ Z:> cr^} 
for a given link obeys Bianchi identities. Integral over IIo-zdo-i d^t'^z would result in 
the singularity of the type of [5^(i?o-2 — -^^2)]^ for some D a^. To avoid this, we 
should omit integration over any one of the triangle containing the link a^. In overall, 
well-defined integral of A/" times monomial of all ^"28 is that over area tensors of the 
triangles not contained in some set J^, no-a^j^d^Vo-a. It is not difficult to see that this 
is the same set T of the triangles considered in the Introduction on which curvatures 
could be expressed in terms of the curvatures on other triangles. Generally it can be 
defined as the set of the triangles with edges constituting a set of non-intersecting and 
non-sclf-intersecting unclosed broken lines (that is, only two such edges can meet at 
any vertex) passing through all the vertices of our manifold. We can refer for a moment 
to this set of finks as to the discrete analogs of the lapse-shift vectors which pick up 
direction of a coordinate t at each vertex. So the triangles might be called "i-fike" 
ones. Omitting integration over some area tensors is equivalent to saying that these 
are fixed. This resembles the issue of the gauge fixing in the continuum theory. There 
is even direct relation to such gauge fixing. Indeed, we can pass to the continuous time 
theory by shrinking the edges along any direction temporarily chosen as time. In this 
limit we can develop Hamiltonian formalism and find [18] that fixing area tensors of the 
triangles of is an admissible gauge fixing in the resulting theory. This is analogous to 
fixing lapse-shift vectors in the continuum general relativity with additional difference 
that now when studying analytical properties of J\f, as mentioned in Introduction, area 
tensors are treated as formally independent variables generally not corresponding to 
certain lapse-shift vectors, therefore fixed are namely area tensors of the i-like triangles. 
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not lapse-shift vectors. 

Thus, the moments of the distribution M can be defined as integrals of M with area 
tensor monomials over d^v^r'^. ^ ^ T . We have seen this for a simplified form of action 
with g{x) — X. At g{x) — arcsin x restrictions on the set of area tensor components 
to integrate over may be relaxed (since instead of ^^{Ra-i — -RJ2) we have less singular 
distribution when integrating exponent over d^fo.2), and this set might be larger than 
{v„2 : (T^ ^ JF}; this point requires further studying. In any case, omitting integrations 
over d^t'o.2, cr^ e ^ is sufficient to avoid singularity when integrating exponent over 
^Vfji%. Next, it is more illustrative to redenote area tensors of the triangles of ^ as 
v„i — > T„i to stress that these are temporarily referred to as constants. Contribution 
of the triangles of T in exp( jS') decouples multiplicatively as 



1 + - I \/ "^T^2 arcsin 



V,2 * +i?,2({i?,2 : a2 ^ ^}, {f],3 : G Q}) 



+ 1 \/ '''^2 arcsin 



-T,2*-i?,2({i?,2 : a2 ^ ^},{Q,3 : e Q}) 



n^^-^ (40) 



where i?o-2, cP' e T ^ is written by Bianchi identities as function of other -Rct2s, ^ T 
and ilo-3S, G Q. The earlier neglected dots in (19) are just due to (40). This is 
analytical function of =^^2 = e^bc^^^V^, analytical at ='r^2 = 0, cr^ ^ J^. 

Finiteness of the integral of a function with any product of its arguments means 
that this function is decreasing faster than any inverse power of arguments. The 
simplest such function is exponentially decreasing one. This type of decreasing is 
most natural in the present case when the function of interest H is itself integral 
of exponent (and, by proper deformation of integration contours in complex plane, 
is expected to be representable as a priori combination of increasing and decreasing 
monotonic exponents). Exponential suppression should take place at least over area 
tensors of a part of triangles when area tensors of other triangles T (for which we 
can take t-like ones for any coordinate t) arc fixed. In the particular case r^'i = 0, 
if (7^ e J^, we simply get factorization of Af into the above calculated exponentially 
decaying factors for separate triangles cr^ ^ J^. 



4.4 Suppression of the edge lengths 

Examining the measure over link lengths we find suppression factors for degenerate 
(zero volume) configurations of spacetime, and, on the other hand, the number of 
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triangles suppressed at large areas at present dimensionality four of spacetime turns 
out to be large enough to provide also suppression of large edge lengths. 

In more detail, consider area tensor part of the path integral measure. The sets of 
variables {v^^ : cr^ ^ J^} and {Jl^s ■ cr^ ^ G} become true canonically conjugate sets in 
the limit when some coordinate (not necessarily t) is made continuous (by shrinking 
simplicial sizes along it) and regarded as time in the Hamiltonian analysis [18]. We are 
trying to define the full discrete measure which would become true canonical one of the 
type of dpdq in the continuous time limit. As usual for discrete counterpart of a con- 
tinuum theory, some first class (commuting) constraints convert into non-commuting 
second class ones which, however, now remain first class in the fiat spacetime. Then 
the constructed in the standard way canonical measure is singular at the fiat space- 
time. The way round this is to construct the measure in the theory with independent 
area tensors which, accidentally, have proved to be good also for defining integrals 
over connections in the present paper. As a result, we get a theory with first class 
constraints (at least, in empty spacetime). The constructed measure is simple and has 
the area tensor part Yia'^^j^d^'Vo-^- Then the measure could be projected onto the hy- 
persurfacc of genuine gravity in the configuration superspace singled out by geometric 
relations between area tensors required for that the edge vectors would exist on which 
area tensors could be constructed as bivectors. 

To illustrate this in the notations somewhat close to the continuum ones using 
tetrad, we again assume certain ordering of the simplicial structure w.r.t. a coordinate 
t direction, see fig.l. Then each 4-simplex has a t-like edge (i. e. of the type of (i^i))- 
A 4-simplex is spanned by a tetrad in the 4-simplex local frame emanating from some 
of its vertices. Let 1^ be 4-vector of the t-like edge, an analog of the lapse-shift vector. 
The bivectors are v^^^ = ^^"■^ cdl^lfi- The triangles having v^^, a = 1, 2, 3 and analogous 
in other 4-simplices form the set The area tensors = ^ta^^'^v'^t^, «, /3, 7 = 1, 2, 3 
are those over which integration element 

rid'X'-IldXdX (41) 

a a 

enters the measure. Here d^vd^v* = 2^d^Re'yd^Im v. To confine ourselves by only 
those which are indeed bivectors for certain tetrad, we insert 5-function factor 6^{Va- 
V/j — v* •■y^)(det H^Q-'y/?!!)^. This simply enforces the conditions eabcdv'^v'p' — ensuring 
existence of a triad defining v*^. We imply in the standard way the Shwinger "time 
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gauge" [19] /° = 0, q; = 1, 2, 3 so that the solution is 2vi — il2xh, . . . , perm(l, 2, 3), 

The multiphcr of 5-function is singled out by symmetry consideration that the overall 
5-factor be invariant w.r.t. an arbitrary redefinition of the triad Va Aa^vp leaving 
three-dimensional section in the same hyperplane, i. e. with real A^^ . With this factor, 
the integration element reduces to 

|det||v,-'y^|||3/2J]dV (42) 

(up to integration over 3 rotational degrees of freedom of the triad v* ) . 

Note that the measure is often considered as definable up to some power of — det 
IIS'A/ill- Now g\fj, — Ixlaix in the 4-simplex in certain piecewise-affine coordinates x^. A 
positive power of — dot \\g\f_i\\ = 8(/o)^| det \\Va ■ (at /° = 0, a = 1, 2, 3 assumed) 

means the power of | det Hvq, ■ i'/jH | as well. 

Note also that the factors | det ^v^ ■ V/3||p/^ also arise on three-dimensional faces 
when further reducing (in the most symmetrical way) the area tensor part of the 
measure to the physical hypersurface in the configuration superspace of area tensors. 
The matter is that the number of triangles is considerably larger than the number of 
edges. It is seen for the simplest periodic simplicial decomposition when the spacetime 
is divided into 4-cubes and each 4-cube is divided into 24 4-simplices[13] . Such structure 
contains 15 edges and 50 triangles per vertex. Inside a 4-simplex its 10 edge lengths 
locally can be expressed in terms of 10 triangle areas, but without taking into account 
relations between the different areas the different 4-simplices may not coincide on their 
common faces. When expressing area tensors in terms of edge vectors, superfiuous 
integrations inevitably get transformed into scale factor of positive dimensionality in 
V (once reducing the measure is performed in scale-invariant way, by inserting scale- 
invariant 5-factor). Relations between the different areas are very complicated and 
nonlocal. Formally more simple way to take them into account in the measure is to 
extend the set of variables simultaneously imposing simple conditions. Namely, 
insert for each triangle cr^ a full set of its area tensors defined in the local frames of 
all the 4-simplices v^2\„4 containing this triangle, cr^ D cr^, and viewed as independent 
variables. Each area tensor denoted above as v^^ on which the action depends is 
defined in certain 4-simplex a* = cr^(a^), = "^0-21 0-4(0-2). The action does not depend 
on Vfj2\„4, (T^ 7^ a^{a^). So we can identically insert integrations over L~^d^V(r2\^4, ^ 
(7^ ((7^) in the measure as well, restricting each component v^'i^^ to vary in the region 
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[— L/2, L/2], L oo. After such extension of the set of variables the condition of the 
coincidence of six pairs of edge lengths of the pair of the 4-simplices on their common 
3-face (tetrahedron) is written simply as six bilinear conditions for three area vectors 
defined inside these 4-simplices, and v'^: ■ Vp — v'^ ■ v'^. Take the integration 
element (42) and the analogous one with v'^, |det||v'^ • v^H Ha d^'^^L- Iiisert the 
most symmetrical 5-factor 5^{Va ■ Vp — v'^- v'p){(\ei \\Va ■ Vp\\Y. The multiplier of 5- 
function is singled out by invariance of the overall 5-factor w. r. t. an arbitrary 
redefinition of the tetrahedron edges la Ba^lp, l'^ BoPVp with the same B^^ 
(the viewpoint is undertaken that the primed and unprimed edges are the same in 
the world coordinates, while tangential (affine) metric on the common face in the 4- 
simplices g^p = Z„ • Ip, g'ap = C ' ^'p undergoes discontinuity [20, 21] Ag^p = gap - g'ap)- 
The product of (42) and | det \\v'^ ■ t>^|| |^/^ Ha d'^'^^a transforms into 



(up to integration over 3 rotational degrees of freedom of the triad v'^). According to 
the above said, the relatively large dimensionality in v of this factor is eventually due 
to the approach adopted proceeding through the extended configuration superspace of 
(the large number of the components of) independent area tensors where the latter are 
true canonical variables. The integration over connections of the present paper relies 
on this concept as well, but indirectly, as possibility to analytically continue the result 
to independent area tensors. 

Turning to the link suppression, there is possibility that the triangles with small 
area but the two sides of large length (spikes) might be dominating as is the case in the 
two-dimensional gravity[22]. In our case, as considered in subsection 4.3, exponential 
suppression takes place over v^, a = 1, 2, 3 (area tensors VQa being considered as fixed 
arguments). This means suppression over Zq,s as well, since small v^s mean small Zq,s, 
with exception of the case when |Zi xl2-lz\ = 2\/2\ det ||Va-t^/3|||^^'^ (tetrahedron is 
flattened). The latter case is suppressed by occurrence of the factor | det || • |^/^ in 
the integration element (43). With this factor, the powers of edge vector components 
la — \/2ea^^vp X v.y\ det II Vq-V/jII \~^/'^ have convergent at det || Vq-?;^|| expectations 
< > at n < 20. Were we in two dimensions, we would have two edge vectors of the 
triangle Zi,Z2 and only one area vector li x I2 (more exactly, scalar). Suppression of 
the latter would not be sufficient to guarantee suppression of Zi, ^2- 




(43) 



a 
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4.5 Local maxima 

Appealing feature of the definition of integrals through moments adopted is that when 
computing any moment we only deal with values of finite number of derivatives of g{x) 
at X = 0, i. e. with local properties of 'arcsin' function at x = at each step. This is 
important since thus far path integral formalism has been checked in physics only on 
perturbative level. Nevertheless, knowing full infinite set of moments somehow takes 
into account analytical features of this function. In particular, we observe poles of 
Afoi^: V*) (36) at v^= 4n^(l + i/7)~^, n= 1, 2, .... Corresponding dependence on v looks 
as result of some summation in the path integral over branches of the 'arcsin' function, 
as if we had substituted arcsin arcsin +27rn in the exponential and summed over 
n for each of the two 'arcsin' functions. This would just result in the hyperbolic or 
trigonometric function in the denominator of (36). Thus, knowing full set of moments 
reproduces features related to the global analytical properties of the 'arcsin' function. 

In the physical region Imv^ = proximity of the poles at v^ = 4n^(l + i/7)~^ 
shows up especially at 7 ^ 1 or at 7 ^ 1 as local maxima of A/o(v, v*) approximately 
at v^ — — 47^71^ for 7 1, n -C 7~^ or at v^ — An^ for 7 1, n ^ 7. It is 
interesting to compare location of these maxima with eigenvalues of the operator of v 
(in the continuous time limit formalism). In Ref. [23] we have considered kinetic term 
Tr(vQ'^Q) (3-dimensional version of this one was deduced by Waelbroeck[24] in lattice 
gravity from symmetry considerations) in the canonical form of Regge action. Now we 
only have to rewrite this for nonzero 7"^, 

von-n ^ [1 + -J — \ — ^ + " 7J — 2 — = r " tJ ^^^^ 

for combination of representations +5", ~S (13). Consider the following degrees of 
freedom in Q, 

a— + l3—\ (45) 

(rotation around the triangle parameterized by a and in the triangle parameterized by 
(3 preserving the area tensor). Since v * v = for physical contribution of Q!,/3 to 
the kinetic term is 

For the spacelike area = — |v^|. Then (3 is Euclidean angle, and requiring the 
canonical commutation relations and periodicity at /? ^ + 27r gives eigenvalues of 



20 

the conjugate variable |v|/7: |v| — 771. For the timehke area = |v^|. Then a is 
Euchdean angle, and requiring the canonical commutation relations and periodicity at 
q; — > q; + 27r gives eigenvalues of the conjugate variable |v|: |v| = n. 

It is seen that the v ensuring local maximum in the probability distribution of 
areas turns out to be approximately the eigenvalue of the operator of v within operator 
approach to the continuous time simplicial gravity. This should not seem surprising 
since both effects arc caused by similar reasons, multivalent nature of 'arcsin' function 
and periodicity over (Euclidean) angles. 

5 Conclusion 

To resume, despite exponential growth of the Haar measure, integration over connection 
in the path integral in the discrete simplicial gravity in Minkowsky spacetime can be 
well defined. The result shows that contribution of large areas is suppressed exponen- 
tially. This provides also suppression of large edge lengths which is important for that 
the minisuperspace system described by elementary lengths/areas be self-consistent. 
The set of local maxima in the probability distribution for area has relation to the set 
of eigenvalues of the operator of area in the continuous time theory. 
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